Mean field theory of localization in an elongated burned fuse model 
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We propose a mean field theory for the localization of damage in a quasistatic fuse model on a 
cylinder. Depending on the quenched disorder distribution of the fuse thresholds, we show analyt- 
ically that the system can either stay in a percolation regime up to breakdown, or start at some 
imposed current, to localize starting from the smallest scale (lattice spacing), or instead go to a 
diffuse localization regime where damage starts to concentrate in bands of width scaling as the 
width of the system, but remains diffuse at smaller scales. Depending on the nature of the quenched 
disorder on the fuse thresholds, we derive analytically the phase diagram of the system separating 
these regimes and the current levels for the onset of these possible localizations. We compare these 
predictions to numerical results. 

PACS numbers: 62.20Mk,46.50.+a,46.65.+g, 64.60Cn, 81.40Np,05.40.-a 



I. INTRODUCTION 

To understand breakdown processes in fragile systems 
with elastic interactions between the elements, and dis- 
order in the material properties, fuse networks are often 
studied d 12. Such simplified models correspond to a 
scalar approximation of elasticity, i.e. retain the pres- 
ence of long-range interactions, and such lattice models 
can be conveniently studied numerically, with the possi- 
bility to control a priori the probability distribution func- 
tion (p.d.f.) characterizing the disorder in the rupture 
thresholds [J. Fuse models allow to study the impact 
on breakdown processes, of parameters as the disorder in 
material properties, and of size effects (ratio of system 
size over lattice spacing, or over grain size for a natu- 
ral system). We will present here a detailed study of 
an elongated fuse model, and show how three different 
breakdown regimes are accessible to it depending on the 
nature of the quenched disorder (q.d.) in the rupture 
thresholds, and on the system size. 

Related studies have already been performed on square 
fuse models 0. The present work extends these stud- 
ies to the case of rectangular systems, with an extent 
Ly importantly exceeding the dimension L in the direc- 
tion perpendicular to the main current flow. This ex- 
tension will allow to show how three types of breakdown 
processes can emerge in it, which will be termed as an 
entirely localized regime, a diffuse localization regime, 
and a percolation-like one. We will develop an analytical 
mean field theory, allowing to classify which regime dom- 
inates the final breakdown, as function of the system size, 
and of the characteristics of the quneched disorder. The 
three possible regimes are illustrated in Fig. ^ The total 
localization regime corresponds to breakdown propagat- 
ing between close or nearest neighbors. The percolation- 
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like regime corresponds to systems where a significant 
fraction of the entire set of bonds have to fuse before the 
system becomes nonconducting. The diffuse localization 
regime corresponds to a system where the burned fuse 
concentrate in a band of size comparable with the sys- 
tem width, but where the damage is distributed diffusely 
inside this band, without propagating necessarily to the 
near neighbors of the already burned fuses. 

An important motivation of this study is to charac- 
terize the scaling law between the system size, and the 
characteristic size where damage localizes in the so called 
"diffuse localization" regime. This scaling law has an im- 
portant theoretical impact on the understanding of the 
origin of the geometrical characteristics of natural frac- 
ture surfaces. Indeed, in general, the main contribution 
(so far) to the science of fracture by the physics commu- 
nity over the last twenty years is the discovery that brit- 
tle fracture surfaces are self affine |^ . Self affinity implies 
statistical invariance of fracture surfaces under rescaling 
of length scales parallel to the average fracture plane by 
a factor A and rescaling of the out-of-plane length scale 
by a factor A^, where C is the Hurst or roughness ex- 
ponent. In 1990, based on experimental investigations 
of brittle aluminium fracture surfaces, Bouchaud et al. 
[3 proposed that the roughness exponent has a universal 
value close to 0.8. This value has been reported in maiw 
later investigations, see e.g. Refs. H, l^j 

a small-scale regime governed by a different roughness 
exponent was reported in addition to the "usual" regime 
characterized by a roughness exponent of 0.8, see 
for a review. There have been several attempts at find- 
ing a theoretical explanation for the universal roughness 
exponents, see jlllllllllll- 
Using the fuse model as paradigm for brittle fracture 
P, 12, Hansen and Schmittbuhl jl^] have recently pro- 
posed that the roughness exponent ( is related to the 
exponent v controlling the divergence of the autocorre- 
lation length of the emerging damage, ^, as function of 
the control parameter: more explicitly, in the case of a 
burned fuse model, noting V the imposed global volt- 
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Figure 1: Configuration of burned fuses in an elongated network at system breakdown, for 5 realizations with a decreasing 
quenched disorder from (a) to (e): the distribution of the fuse thresholds t is of the type p{t) ~ t"~^ for < t < 1, with a 
indicate. In the white region underlies a nondisplayed diamond lattice of intact fuses, inclined at 45 degrees with respect to 
the bus bars at the top and bottom, with a lattice step £ = 1. Burned fuses are marked by gray squares, with a gray index 
turning from dark to light in chronological order. In cases (d) and (e), order dominates and the rupture proceeds almost always 
via nearest neighbors: this is an ordered rupture, with a total localization of damage. In case (a), disorder is large and the 
breakdown process is dominated by the distributed location of weakest flaws: this is a percolation-like filling, with a finite 
fraction of bonds to break to reach the system breakdown. In cases (b) and (c), rupture does not proceed via nearest neighbors, 
and looks diffuse at scales below L = 30^, the lateral a;— size of the system. Broken bonds are nonetheless localized in a band 
of vertical size w comparable to the horizontal width L of the system: this is the diffuse localization regime. . 



age difference, and Vc the voltage at complete electrical 
failure , £, ~ \V - Vd'" ■ Hansen and Schmittbuhl [l3| 
proposed the existence of a scaling relationship between 
these two exponents, C = 2i^/(l + 2i^) in such breakdown 
problems. This relationship was numerically checked for 
a fuse model in two dimensions, where both the expo- 
nent v was numerically measured, and found to be the 
one of percolation, v = 4/3 as well as the roughness 
exponent ( — 8/11. The same reasoning for brittle frac- 
ture, based on the same scaling relationship, and an auto 
correlation divergence exponent j/ = 2 [l^, leads to the 
roughness exponent = 4/5 for brittle fracture, in ex- 
cellent agreement with the experimental measurements 
0. Central to this theory is the scaling law w ^ L/£ 
between the width of a concentrated damage zone w and 
the size of the system L, where £ is the lattice constant. 
One of the aims of this paper is to explain the origin of 
this scaling law. 

The fuse model consists of a lattice of ohmic resis- 
tors with identical conductances placed between two bus 
bars, where each bond carries an electrical current up to 
a threshold t above which the bond burns irreversibly. 
Each of these local random thresholds t are fixed ini- 
tially and taken independently of each other from the 
p.d.f. p(t), which entirely characterizes the uncorrelated 



quenched disorder present in this system. 

For square systems, the phase diagram of this system 
was established numerically and through order statistics 
arguments 0] depending of two parameters a and f3 char- 
acterizing the quenched disorder distribution tails in the 
limit of zero or infinite thresholds, as p(t) ~ where 
t — > and p{t) ~ t~f^~^ when t ^ oo. We will con- 
sider here such systems in a cylindrical geometry, i.e. a 
periodic band of finite width L ^ £, where £ is the lat- 
tice constant placed between two bus bar at distances 
Ly ^ L, and derive analytically the equivalent of this 
phase diagram as function of a and L/£, at Ly/£ 3> 1 - 
we will only consider here power-law distributions with 
an upper cutoff, corresponding to /3 — > cx) in the pre- 
vious terminology. With respect to this previous work, 
we extend the study in two ways: we consider elongated 
systems, and we in details the anisotropic aspect of the 
current perturbation generated by burned fuses. 

In the next Section, we present the basic assumptions 
for and philosophy of our statistical analysis of the fuse 
model. In Section Ulll we calculate the shape of the cur- 
rent distribution around a region of burned-out fuses. We 
then present in Section llVI the spatial probability distri- 
bution of subsequent fuse burn-outs. The main result 
of the calculation is presented in Fig. ^ Depending on 
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the disorder exponent a, and on the system size, there 
are three possible breakdown regimes: (1) A percolation- 
like phase where no localization occurs and where a finite 
fraction of the total number of fuses needs to burn out in 
order for the conductance of the lattice to drop to zero in 
the infinite-lattice limit; (2) a diffuse localization phase 
where a damage zone develops, with a width w propor- 
tional to the width of the lattice L; and (3) a complete 
localization phase where a single crack evolves without 
damage around it. These regimes are illustrated in Fig. 
1. We do not in this paper discuss the phase diagram 
with respect to the second disorder exponent /?. In Sec- 
tion we compare our analysis to numerical results on 
the fuse model. We summarize our findings in Section 

ED 



II. MODEL UNDER STUDY AND BASIC 
ASSUMPTIONS 

At any stage of the rupture process, we will assume 
that the local currents in the fuse model are determined 
through a continuous approximation, as the solution of 
the conservation of charge V-j =0 under boundary con- 
ditions j jeH when y ±Ly/2 (the band is L-periodic 
in the a;-direction, are the unit vectors). The cur- 

rent density is of the form j{r) — —c{r)W(j){r), with a 
conductance c(r) equal to unity in the intact cells, and 
zero in the broken ones. 

After the first fuse has burnt at a certain current level 
je at a position defined as the origin, we are interested 
in the average change of external current necessary to 
break the next element: since the problem is linear, for 
a given geometry of burnt elements, the current flow 
for any other value of the external current j' is simply 
(j'/je)ji(r). For a given realization of the quenched dis- 
order t{r) (such as t{r) > je at every location), the next 
fuse will burn when a first threshold is reached by the 
local current, i.e. when the external current reaches 
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at a position r„ corresponding to the realization of this 
minimum. If j„ > je, the applied external current has to 
be increased by a finite value for the next fuse to burn. 
On the contrary, if j„ < jej there is an avalanche and 
the next fuse burns immediately if the external current 
is not reduced immediately during the first burn-out to 
this lower value j„ . 

We are also interested in the geometric characteristics 
of the relative position of the next burnt fuse with re- 
spect to the first one: over all realizations of the quenched 
disorder, we define the probability distribution over this 
relative position of the next burnt fuse as /Lt(r„). Three 
scenarii will be shown to happen, depending on the ran- 
dom mean square distance of the next burnt fuse to the 
previous one, cP = J r^^(r„)dr„: (1) d ^ +oo and the 
process remains diffuse, resembling a percolation process. 



(2) d ^ £ i.e. it is a function of the lattice spacing, in- 
dependent of the system width L. This is the onset of a 
complete localization, i.e. the current perturbation cre- 
ated by the broken cell is such that the rupture will prop- 
agate mainly from nearest neighbor to nearest neighbor 
up to complete breakdown of the system. (3) d L, 
which is the onset of a regime which we define as "diffuse 
localization:" damage starts to concentrate in a band of 
a width in the y direction comparable to the system size 
in the x direction, L, but the closest neighbors of the 
previously burnt cell are not significantly favorized. This 
is the regime where the scaling arguments of 0j should 
apply. 

If the system remains in the diffuse regime, the spatial 
correlations of the damage are not significant, and we 
are entitled to consider a mean field approximation to 
study the subsequent history of the process: if the last 
fuse has burnt at a location at a current level je, the 
probability distribution over the location of the next fuse 
burning is approximated as the probability obtained from 
a situation where a single fuse has burnt at rp, under 
the condition that all of the remaining thresholds where 
above je ■ 

To estimate the average level of current necessary to 
trigger the next fuse burning and the statistical proper- 
ties of its location, we extend the arguments of Roux and 
Hansen {TjJ: By convention, any level of local current j 
in the system will be expressed through a reduced di- 
mensionless variable s = (j — je)/ je, ratio of the current 
perturbation generated by the last fuse burnt, over the 
average imposed current level. We next define n(s)As as 
the number of cells experiencing a local current between 
je(l + s) and je(l + s + As), where As <C 1 is a small 
parameter. Defining as f2(s. As) the region experiencing 
that local current level, we have 



n(s) — lim 



1 



As^O ^2. As 



dxdy 



{x,y)eft{s,As) 



(2) 



The average value m of the external current leading to 
the next burn-out is, from Eq. the average value of 
the minimum over all cells of the random variable y = 
t/{l + s) — Eq. O, 



^inin{s=nAs,ren(s,As)/«6Z} ^ 



t(r) 



(3) 



At a given level of current perturbation s, we define 
P{y, s) as the cumulative probability of the random vari- 
able y = t/{l + s), given that t > je- This last condition 
refiects the fact that the intact fuses experiencing a cur- 
rent je{l + s) have survived, up to the burning point 
of the fuse creating the dipolar perturbation we look at. 
This is straightforwardly 

Piy,s)^ ^ _ p^. ^ Re{y{l + s) ~ Je) , (4) 

where He is the Heaviside function, and P is the cumu- 
lative distribution of thresholds. 
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As shown in the Appendix. |XI where we extended some 
statistical resuhs of Gumbel [Ig , m satisfies the implicit 
equation 

/ n{s)P{m, s)ds = 1 . (5) 

J s 

We also show in this Appendix that fj.{rn) a P{m, s(r„)), 
where m is the solution of the above equation, i.e. that 
n{s)P{'m,s)As is the probability that the next bond 
would break in il(s, As). Thus, if we find {smax, As) such 
as the integral in the above has a significant support only 
in [smax, Smax + As] - i.e. /j^^T n{s)P{m, s)ds = 1., 
the next break will almost certainly happen in the spa- 
tial region Q{smax, As), and the geometric properties of 
this spatial ensemble are representative of the ones of the 
spatial distribution over all possible locations of the next 
broken bond, i.e. the random mean square distance to 
the next broken bond will be evaluated as 

^ [ {x^ + y^)dxdy. (6) 

III. NUMBER DENSITY OF CELLS OVER THE 
LEVEL OF CURRENT PERTURBATION. 
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Figure 2: Contour map of the elementary current perturba- 
tion due to a burnt fuse at the origin. The long-dashed curve 
corresponds to a 0-perturbation, the short-dashed one meets 
a saddle point at (7r,0), and corresponds to f{u,v) — 0.5. 



We wiU now compute the mass n(s)As and shape 
i7(s. As) of each region carrying a certain value of the 
local current magnitude in [je(l + s),je(l + s + As)]. 
The local current, after an unit has fused somewhere, 
is written as j{r) — j^y + Sj{r), with a perturbation 
Sjir) = —Vcj) and a potential field satisfying Laplace 
equation = under Neumann boundary conditions, 
V0 = when y — )■ ±Lj, and fi ■ \J4> = Jen • y along the 
surface of the broken element (elementary lattice cell), 
where fi is the elementary vector normal to it. Since 
L Ly, this current perturbation will be approximated 
as the one in an infinitely long cylinder, i.e. the long 
range condition used will be V0 = when y — > ±c», and 
x-periodicity with a period L. We will then also use 
the coordinate system where the last burnt fuse ise at 
the origin. Furthermore, from a distance of a few lattice 
size and above, the shape of the lattice cell is no more 
relevant, and this elementary current perturbation is it- 
self approximated as the solution of this problem with 
a spherical fused element of diameter i: (p satisfies in 
circular coordinates, h-V(j){r = i/2,9) — jeSm{d). For 
sufficiently large systems L/£^ 1, this particular poten- 
tial can itself be constructed as 



-je7re^y-VG/2 



(7) 



where G is the solution of the Poisson equation in 
periodic Boundary Conditions in the x direction, sat- 
isfying AG = 5(x,y) and G(x + L,y) — G(x,y): in- 
deed, along the surface of the elementary circle of di- 
ameter we have VG ~ f/27rr, and with f the ele- 
mentary radial vector, and the angle between x and f. 



f-\I{y-WG) = f-V(sin(6l)/27rr) = -2 sin(6')/7r£2. The 
complete expression of G in such periodic boundary con- 
ditions is after Morse and Feshbach [191, 



4sin'^ 



!I^U4sinh^ 



L I 



V L 



(8) 



Eventually, at a sufficient distance fro m a broken c ell 
r 3> ^, we have 5j ^ je and j{r) ~ yj (j^y -t- 5j{r) ~ 
je + where 5i{r) = y-5j{r) - which is a classical 

expression for the dipolar perturbation emanating from a 
burned fuse in such models, see e.g. |23|. The magnitude 
of the current perturbation is thus determined from the 
above Eqs. (|7I8(I as 



Sj{x,y)/je 
with f{u,v) 



^f{27:x/L,2ny/L) 

1 — cos (u) cosh (v) 
(cosh (v) — cos (u))^ 



(9) 
(10) 



A contour map of the dimensionless current perturbation 
f{u, v) is displayed in Figure |2 

Since this perturbation /(u, v) is pair in both its argu- 
ments, only a zone [0, tt]^ was represented. The system is 
27r-periodic in the x-direction. Two special contours were 
highlightened: f{u,v) = is the long-dashed curve. On 
the displayed region, points to the right of this line expe- 
rience an increased current due to the burnt fuse at the 
origin, and conversely the current is screened for those to 
the left of it. This zero perturbation contour correspond 
to Vo{u) = acosh(l/ cos(u)), which has a support on 
u mod [27r] G [— 7r/2, tt/2] and an asymptot vo{u) — > +oo 
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Figure 3: Dimensionless number density of cells as function 
of the level of current perturbation, and asymptotic forms in 
dashed, at infinite distance (a=0), around the saddle point 
(a=l/2) and in the region close to the burnt fuse {a oo). 

when u — > ±tt/2. The other contour goes through a 
saddle point of / in (u,v) — (0,7r), and corresponds to 
f{u,v) = 0.5, or vo.5{u) = acosh (\/2 — cos^ u ) . 

The regions fl{s, As) that we want to characterize geo- 
metrically, which support a perturbation of current such 
as s < 6j / je < 5 + As, correspond to the regions between 
two neighboring lines of the contour map on Fig. |21 The 
number of the cells in such regions, defined in Eq. ||2Jl, is 
shown in Appendix IbI to be of the form 

'^(^) = ^9i2Lh/nH^), (11) 

where the dimensionless quantity g, function of its dimen- 
sionless argument, is numerically evaluated and plotted 
as a continuous line in Fig. |3 The numerical evaluation 
is based on analytical expressions detailed in Appendix 
IbI where are also derived the asymptotic behaviors of 
this function g which are plotted in Fig. |3 

IV. REGION OF MOST PROBABLE NEXT 
EVENT 

We have now derived the number of cells associated 
with each current level, n(s), characterizing the interac- 
tions in this system, and need to consider some specific 
quenched disorder to determine the typical separation 
between two subsequent burning fuses. With a quenched 
disorder distribution of power law type P(t) = on 
Q < t < 1, with a > 0, we obtain the cumulative distri- 
bution for thresholds to be below m, for the fuses that 
were still intact at current je, through Eq. as 

P{m,s)= ^ ;'[^ ^- He(m(l-Ks)-je). (12) 

We divide the space with respect to the last burned 
fuse in three zones, noted f2c, and D,f, and have to 



solve the implicit Eq. jSJ to find both the most proba- 
ble region out of these three where the next break will 
happen, and the most probable value of the external cur- 
rent m at which the next fuse will burn. This equation 
becomes 

(/ +/ +/ )n{s)P{m,s)ds = l (13) 

By definition, ilc is a region of fuses close to the last 
burned one, with the largest positive current perturba- 
tion, such as /(u, v) > Lj (27r^^). All such fuses lie within 
a distance r < \f(L£) of the last burned one, where we 
recall that i is the lattice step, and L is the a;— dimension 
of the system. 

r2(i is a region defined with moderate and finite current 
perturbations, where Lji^TV^l) > f{u,v) > 1/4. The 
typical distance r from the origin of the current pertur- 
bation, over the zone Qd, is such as y^L£) < r < ttL. 

Last, is a region of weak to negative current pertur- 
bation, defined by the implicit equation f{u, v) < 1/4. It 
will be shown that when this region dominates the left 
hand side of Eq. H13|l , the leading contribution to it comes 
from points sitting at a characteristic distance r to the 
last burned fuse, scaling with the system size as r Ly. 

In Appendix El we analyze in details the three terms 
of Eq. H13|l , and reformulate it as 

iJe(A) + Ha{\) + Hf{\) = (14) 

Je 

where A — m/je, with jej"^ the values of te external 
current at the last break and at the most probable next 
one, £ = 'n'^ t"^ / {2L?) and Hd, Hf are proportional to 
the integrals in Eq. (|13|l over the regions ilc, ^d, ^f- 

We will classify the regime of the system according to 
the dominant term in the left hand side of Eq. (|14|l : If 
Hf dominates, the system remains in a diffuse regime 
where there are no noticeable spatial correlations in the 
pattern of burnt fuse. If He dominates, this signifys the 
onset of a complete localization regime where the dam- 
age will develop in a concentrated zone scaling as the 
lattice size £, and propagate through the system, tear- 
ing it with jumps between successive events close to this 
smallest scale. Last, the dominance of Hd would denote 
the onset of a diffuse localization regime, where the char- 
acteristic distance d between the burnt fuse scales as L, 
the system's width. 

Thus, following as the imposed current increases, 
which of the three terms dominates in Eq. (|14|) , allows to 
understand when damage starts to localize, and at which 
spatial scale. This allows to classify, as function of the 
system dimensions Lx/(, Ly/£ and of the quenched dis- 
order, characterized by a, in which localization regime 
the system ends up in. 

It is shown in Appendix [C] that in the early stages 
of the process, at small je, Hf dominates the solution 
of Eq. H14|) . owing to the singularity of n(s) around 
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s ~ (zero current perturbation line), and this equation 
reduces to 



A" - 1 



1 



Ncellsje 



(15) 



Since the first break is typically for j" — 1 /N cells , this 
equation predicts a second break typically at 



A? -1-1, 



(16) 



i.e. the second break should happen on average at j2 = 
Aiji = (1 + l)i/"ji = {2/Nceiuy/''. Since j2 > Ji, 
the process is stable and there is a finite gap in external 
current to trigger the next fuse burning. Since Hj is 
dominated by the asymptot of zero current perturbation 
(noted hi and he in Appendix C), corresponding to the 
long dashed curve in Fig. [3 which spans the whole y— 
range of the system, the next fuse is likely to burn at a 
distance scaling as d ^ Ly from the first one, i.e. the 
system remains in a diffuse regime, with no noticeable 
correlations between the locations of the burnt fuses: the 
size of the system wins compared to the attractive feature 
of the current concentration around the last burnt fuse, 
in a Flory-type argument. We can then proceed with this 
mean-field theory to treat the later stages of the process. 

As long as Hf dominates Eq. H14|l . Ea. H15|) remains 
valid, and by recurrence, we show in AppendixlOtha-t the 
n — th fuse burns on average when the external current 
is such as 



j: = (n/NceUs) : 



(17) 



As long as this is the case, the n-th weakest bonds are 
the most likely to be the n first burned ones. 

For threshold distributions characterized by a very 
large disorder, i.e. in the limit a —> 0, it is shown in Ap- 
pendixElthat Hf{X) always dominates in the solution of 
Eq. H14() . up to the moment where = 1/2. In this limit, 
the n— th fuse burning corresponds to the n — th weakest 
threshold, and this lasts until the entire system is broken 
due to burned fuses percolating through the system. In 
this case, the process remains diffuse, in a percolation- 
like regime, up to the moment where P{je) = 1/2, which 
corresponds to the critical percolation threshold. This 
means that in this limit of non-renormalizability of the 
q.d. distribution, and very large disorder, the process 
is equivalent to a bond-percolation process, which was 
shown by Roux et al. 21]. 

On contrary, for very small disorder, in the limit a — > 
-hoo, we show in Appendix lUl that Hc{\) dominates, and 
even that the contribution of the nearest neighboring cells 
on the sides of the last burned fuse dominate the integral: 
thus, the next fuses to burn are the ones carrying the 
highest current perturbation, and from the asymptotic 
expression of Hc{X) derived in AppendixlUl and Eq. (|14|l . 



the level of next break is set by 

2l/Q 



A 



where s{a) 



jes(a)i/" 



272 



c?7. 



(18) 
(19) 



This happens in a controlled way, i.e. for A > 1 if a 
is still sufficiently small so that s{a) < 2N — 2/j" , or 
through immediate avalanches A < 1 in the opposite case. 
This is the limit of no disorder, where all bonds share the 
same threshold, and the concentration of current around 
the first broken one is the significant parameter control- 
ling the process in this case: the rupture proceeds from 
the smallest scales, expanding through nearest neighbors 
from the initial seed to tear the system apart. This corre- 
sponds to a classical rupture process, analog to the rup- 
ture of a perfectly elastic and homogeneous material (no 
disorder), where the stress concentration at the tips of 
an initial default leads to the rupture of the system when 
the load is increased - in stable way or not, depending on 
the load level - the situation known from one century in 
linear elastic fracture mechanics, treated by Griffith and 
Inglis 13. 

Between these two extreme cases, in the range of fi- 
nite a, the system can be driven to a third regime if Hd 
dominate in the solution of Eq. (|14|l : correlations in the 
damage start to be significant, but the characteristic dis- 
tance to the preceding burnt fuses is in a range between 
^/L£ and L, and does not scale as the lattice constant i: 
this is the regime which we refer to as "diffuse localiza- 
tion" . 

We determine a lower value am of the exponent of the 
q.d. distribution separating systems entirely equivalent 
to percolation up to breakdown, and these leading to 
diffuse localization, as follows: as long as the percolation 
regime holds, the value of external current, and the size 
of the jumps A in it, are determined by Eq. H17|) . This 
regime goes on as long as Hd{X) can indeed be neglected 
in front of Hf (A). If both terms become equal, the system 
transits towards the diffuse localization regime, which is 
shown in Appendix |C1 to correspond to leading order in 
l/Nceiisi to the condition 



^ln(^:^ 

2 ^e' 



1 - 



(20) 



If this condition is not met at percolation threshold 
1/2, i.e. if 



a < a„ 



ln(L/£)' 



(21) 



the system always remains in the percolation universal- 
ity class. If on contrary a > a^n, the system undergoes a 
transition towards diffuse localization at a typical exter- 
nal current 



J, = l/(H-aln(L/^)/4)i/«. 



(22) 
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Similarly, we determine an upper cutoff aM of the ex- 
ponent of the q.d. distribution, above which complete lo- 
calization will prevail about the diffuse one. By equating 
Hf{X) and Hc{X), with A evaluated from the percolation 
regime expression in Eq. H15() , we show in Appendix El 
that this upper cutoff satisfies the implicit equation 

sjaM) - 27rL/l ^ \n{L/() ^^^^ 
aM 2 

From the expression of s, the equation above has a single 
solution at aM = 1- If a > 1, the system will transit 
towards complete localization at a characteristic current 
level 

jt - 1/(1 + {s{a) - 2^L/^)/2)i/". (24) 

Eventually, ii L/£ < e'^ ~ 54, the above would lead to 
aM < dim 1 and no diffuse localization is obtained for any 
value of a. There is instead a transition directly from 
the percolation regime below a < to a regime leading 
to complete localization for a > a^, with 

o-TT r TT 1 

sK) - ^ = 7 ln(— )(A^_o - 1) = 2. (25) 

where Xperco is estimated by Eq. p5|l . For a > ad, the 
system starts a complete localization at a current level 
given by Eq. . 

To summarize the above results, a phase diagram of 
the system, showing the regime through which it will go 
to final breakdown, is shown in Figure 2| The value of 
ad was determined numerically from Eq. 1)25(1 . A visual 
representation of sequences of burning fuses, for small 
systems, in 5 points of this phase diagram, illustrating 
the three regimes, is given in Fig. ^ 

The above can be compared, in the limit of in- 
finitely large systems, to the numerical analysis carried 
by Hansen et al. using the notations of this paper, 
a = 4>o, and l/</>oo = 0, and the system goes from a dis- 
orderless regime A when a > aM = 1 to a scaling regime 
B with diffuse damage and localization when a < a]\j. 
The difference between the critical exponent separating 
the two regimes, which is aM = 1 in the present case, 
and (/io = 2 in the models of |j3| , is believed to come from 
the elongated character of the systems considered here 

{Ly > L). 

Some remarks can be done on the succession of approx- 
imations carried to establish this phase diagram: Most of 
these approximations correspond to keeping the leading 
order in the inverse of the number of lattice cells in the 
lateral dimension, i/L, for systems of infinite anisotropic 
ratio (such as L/Ly —^ 0). These approximations al- 
ready take into account finite size effects, since they keep 
finite £/L. They should thus be valid as long as these 
numbers are i <^ L, and L ^ Ly. Evaluating the follow- 
ing terms corresponding to higher orders of the parame- 
ters (./ L and L/Ly , in order to estimate the quality of 
the asymptotic expansions reduced to leading order, is 




2 4 6 8 10 12 14 
In(Ul) 

Figure 4: Phase diagram of the system displaying the regime 
through which it will go to macroscopic breakdown, as func- 
tion of the system's width and exponent characterizing the 
fuse threshold distribution. 



beyond the scope of the present work. In this analyti- 
cal development, there is however an approximation that 
does not fall in this category of asymptotic expansion: 
in order to evaluate whether the system departs from 
the percolation like regime, as more and more fuses are 
burned and the system has stayed so far in this regime, 
we have considered separately the current perturbation 
triggered by each burned fuse. This is perfectly justified 
in the early stages of the process, where since the pro- 
cess is in a percolation like regime, successively burned 
fuses are at distances of order Ly from each other, and al- 
most do not interact. However, as the density of burned 
fuses increases, it can become finite (for sizes and q.d. 
where the process always remains in a percolation like 
regime, the density reaches eventually a large fraction, in 
principle the percolation threshold). In this situation of 
high fraction of burned fuses, the approximation corre- 
sponding to evaluate the local current as a homogeneous 
background je, superimposed to the perturbation ema- 
nating from a single burned fuse (the last burned one), 
becomes of lower quality, due to the existence of multiple 
close burned fuses anywhere in the system. Overcoming 
this limitation at high fractions of burned fuses, would 
require to take into account a large number of perturba- 
tion sources. This task seems more suitable for a purely 
numerical approach: a main scope of this paper is to 
carry out the analytical development in terms of extreme 
statistics for a mean field theory going beyond the purely 
homogeneous description of damage (i.e. incorporating a 
homogeneous term plus a local perturbation). Carrying 
out analytically the details of the calculation while keep- 
ing the track of every local perturbation is beyond the 
scope of this work. The result of this approximation is 
mainly to underestimate the number of close perturba- 
tion sources anywhere in the system as the process goes 
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on: thus, this overestimates the weight of the probability 
to break far from the existing sources, Hf, and under- 
estimates the probabihty to enter a diffuse locahzation 
or a total localization regime. Thus, the mean field the- 
ory presented here should predict properly the transi- 
tion between "total" and "diffuse" localization regime, 
but should over estimate the domain of the "percolation- 
like" regime: in Fig. the left line should be located at 
smaller sizes. This approximation seems to overestimate 
the transition size L/ihy a finite factor not exceeding an 
order of magnitude, as will be shown in the next section. 

Eventually, we note that in the diffuse localization 
regime, the process looks uncorrelated at the lattice con- 
stant scale, i.e. looks like a percolation system, but the 
arguments developed in this paper show that damage 
starts to concentrate in a band of width scaling as the 
width of the system L. An argument based on percola- 
tion in a gradient corresponding to the structure of the 
damage concentration at the scale of the system can then 
be applied to describe the breakdown process, which sus- 
tains the arguments developed in [l4j to explain the ori- 
gin of the roughness of the ultimate breakdown connected 
fronts in this regime. 

Qualitatively, the phase diagram shown above sup- 
ports the idea that the failure of natural macroscopic 
heterogeneous systems is dominated by either the "to- 
tal localization", or the "diffuse localization" regime. 
Indeed, macroscopic materials are often systems much 
larger that the typical scale of the disorder, i.e. sys- 
tems with a high ratio of system size over cell size, L/ell. 
In such regime, the present work predicts that the per- 
colation regime vanishes. More precisely, in this limit, 
the percolation regime would only subsist in the limit of 
non-normalizable threshold distribution, corresponding 
to q; — *■ 0. So the present work predicts that the break- 
down of such system is "totally localizing" at low disor- 
der, or "diffusely localizating" at larger one. This picture 
is consistent with the fracture properties of natural ob- 
jects: when the fracturing solid is more homogeneous, 
or has only moderately disordered toughness properties, 
corresponding to large values of a, the rupture is initi- 
ated on the weakest flaws, and fracture propagates from 
nearest neighbor to nearest neighbor: this is the classi- 
cal picture of linear elastic fracture mechanics of a ho- 
mogeneous solid, described here as "total localization". 
The fracture of such regular object, as e.g. a crystal, 
leaves a flat, or close to flat, fracture surface, as seen in 
Fig. 1(e). Conversely, when the toughness properties 
of the breaking solid are more scattered, i.e. at smaller 
a, when the heterogeneous solid is more disordered, the 
rupture proceeds according to the "diffuse localization" 
regime (illustrated in Fig. 1(b)): this corresponds to the 
rough post mortem fracture surfaces observed in most 
natural materials, found to be self afhne with a universal 
roughness of 0.8. 



a=l5 




break 



Figure 5: Total localization: distance between successive 
burning events, as function of their index, averaged over 50 
realizations. For a quenched disorder exponent a — 1.5, this 
distance remains of order unity, irrespectively of the dimen- 
sions L and L„ of the lattice. 



V. COMPARISON TO NUMERICAL 
SIMULATIONS 

We now turn to confront this theory to numerical sim- 
ulations of the fuse model. We consider rectangular mod- 
els oi L X Ly cells, with high aspect ratios Ly/L in 
order to be close to the infinitely long cylinder consid- 
ered so far. The lattice constant t is now considered as 
unit length, the models considered are periodic along the 
transverse a;-direction, and the rows of nodes at both lat- 
tice boundaries are set to two constant potential values, 
with a voltage drop AJ7 between both regularly increased 
from 0. The rows of fuses are inclined at 45 degrees 
with respect to the x and y directions. Current conser- 
vation (Kirchhoff equation) is required at each node, and 
the current through each fuse at location r connecting 
neighboring nodes with a local potential drop ISV be- 
tween them, is ISV if the fuse is intact, or if the fuse 
is burned. This allows for each configuration of burned 
and intact fuses, to obtain by solving a linear system the 
voltage at each node, and the corresponding current map 
as j (r)=C(r) Af7, where C(r) depends on the configura- 
tion of burned and intact fuses. The linear inversion is 
performed via a conjugate gradient algorithm (Hestenes- 
Stiefel, Eqs. (32) to (38) in Ref. 23]). Initially, the sys- 
tem is entirely intact and random thresholds of maximum 
sustainable current jt('') are picked from the quenched 
disorder distribution, independently for each fuse. At 
any stage of the process, the location r of the next fuse to 
burn and the corresponding value of the external current 
Af7 is obtained as AC/ = min^ [jt(s)/C(s)] = jt{r)/C{r). 

An example of configurations and history of burning 
fuses is displayed in Fig^ for systems of size L x Ly = 
30 X 100, and values of a between 0.25 and 2.5. The char- 
acteristic features of the three regimes are examplified in 
these cases. 

For various values of L, Ly and of the exponent a char- 
acterizing the quenched disorder, we look at the distance 
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Figure 6: Fuse threshold distribution corresponding to ex- 
ponent a = 0.25: percolation universality class. Distance 
between successive burning fuses, divided by the lattice elon- 
gation Ly, averaged over 50 realizations and 200 successive 
events. Irrespectively of the lattice dimensions, this distance 
is comparable to Ly (of order Ly/ 10 here). Note that this 
corresponds depending on the lattice dimensions, to an aver- 
age distance d equal to 20 to 200 lattice units, as the insert 
shows, and equal to 1 to 10 times L. 



d between two successive events, as function of its oc- 
curence number in the succession of events up to total 
failure of the system (when a connected line separated 
the upper and lower boundaries of the system). This 
distance is averaged over 50 realizations. 

The characteristic situation corresponding to a > 1 
is illustrated for the case of a = 1.5 in Figure [SJ the 
distance between successive events is from a very early 
stage of order of a few unities, irrespectively of the sizes 
L (4, 10, 20) and Ly (200 and 2000) considered. This 
regime was referred above as "total localization" . 

On contrary, for low exponents a corresponding to 
larger disorder, the typical situation is illustrated on 
Figure El by the case of a = 0.25: Simulations have 
been performed using lattice elongations Ly = 200 and 
Ly = 2000, and widths i = 4, 10 and 20. The distance 
between successive events has been averaged over 50 sim- 
ulations. Even so, this quantity is still highly fluctuating, 
and an additional running average over 200 successive 
events is performed in order to extract the proper slowly 
varying average of this distance. This high fluctuations 
are easily explainable: this regime is expected to be in 
the universality class of percolation, where the distribu- 
tion of this distance at any stage is non negligeable for all 
possible distances in the lattice. Assuming consequently 
that the root of the variance of this distribution is of 
the same order as its average, the central limit theorem 
ensures that the root of the variance of the averaged dis- 
tance over 50 realizations is still of the order 1/7 of its 
average, which still corresponds to a high noise to signal 
ratio. The resulting average distance is plotted in figure 
scaled by the lattice elongation Ly in the main figure, 
or directly in lattice constant units in the insert. This 
distance is found out to vary slightly during the process, 
and shows 50% variations between the different probed 
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Figure 7: a — 0.5: features of localization, or of percola- 
tion, depending on L and Ly-. Localization (short distance 
between successive events) is seen for larger L and smaller 
Ly, whereas episodes with distances significantly larger than 
L are observed for short L and large Ly. Data averaged over 
50 simulations. 



widths L, but the main result is that the average d is of 
order O.lLy, i.e. scales with Ly: this is consistent with 
the prediction of the previous sections, that systems of 
infinite elongation Ly are isomorphic to percolation, i.e. 
that the distribution of burnt-out fuses is homogeneous, 
irrespectively of the configuration of the already burned 
fuses - which would predict for a very elongated system 
— !■ cx), an average distance between successive events 

lo" lo" dyidy2\yi - y2\/Ll = Ly/3. The fact that 
we observe d ~ Ly/ 10 rather than Ly/3 can be under- 
stood as a finite size effect: less cells far away from the 
last fuse burning are likely to present the minimum ra- 
tio t/ j, which increases the likelihood of having a next 
burned fuse in the zone of significant current perturba- 
tion, closer to the last burned fuse. 

We have also analyzed the behavior of the system for 
a = 0.5, where according to Eq. ( I21|l . in the limit 
of Lj, — > oo one expects a percolation-like behavior for 
L < e^/°-^ - 2980, or a diffuse localization behav- 
ior for larger system width. The average distance for 
L = 4, 10, 20 and Ly = 200 and 2000 is displayed on fig- 
ure 0- the insert represents the same data on a smaller 
scale. A qualitative interpretation of these results can be 
presented as follows: focusing first on the least elongated 
systems {Ly = 200), the average distance is for L = 10 
and 20 of a few units, but the thinest systems, L = 4, 
displays a more complicated behavior: after an initial 
decrease, the distance displays episodes where its mag- 
nitude is around a few unities, alternated with episodes 
of order Ly. This can be interpreted as a case on the 
verge between localizing or non-localizing, i.e. as a case 
sitting around the line separating percolation from lo- 
calizing regime in fig. ^ The value L = 4 is consider- 
ably smaller than the predicted L ~ 2980 for infinitely 
elongated system. This presumably results from the un- 
derestimate in the analytical calculations, of the localiz- 
ing / nonlocalizing separation due to the high fraction 
of burned fuse close to the breakdown process, and from 
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Figure 8; Distinction between total and diffuse localization: 
at elongations Ly = 200 or 2000, average distance d between 
successive events, for the second half of events, as function 
of the lateral size L of the system, on a bilogarithmic scale. 
Over the numerically accessible range as L grows, d is seen to 
saturate for a = 1.5, corresponding to total localization. For 
a = 0.75, djL^I'^ does not vanish , corresponding to diffuse 
localization. 

strong finite-size effects at finite Ly^ as explained in the 
previous Section. 

The presence of important finite-size effects is con- 
firmed by the fact that for more elongated systems, 
Ly — 2000, such episodes where the average distance 
significantly exceeds the width of the system occur even 
more for the case i = 4, and appear also for the case 
L = 10, while they are absent from the case L = 20: 
presumably, the boundary between localizing and per- 
colating systems is around L = 10 in that case. This 
means at finite elongations Lj,, this boundary for any 
a corresponds to significantly smaller L than value pre- 
dicted for infinitely elongated systems. This finite size 
effect is observed to diminish for increasing elongations, 
as expected: the larger is the elongation Ly, the larger 
is the transition width L. Due to numerical costs, it 
seems however difficult to evaluate numerically, in the 
limit Lj, — !■ oo where this finite size effects would vanish, 
the exact transition value for the x-size separating non 
localizing systems and systems with diffuse localization, 
for example for such q.d. at a = 0.5. From the above, at 
finite Ly, this transition size &i a — 0.5 is bounded be- 
tween L = 10 (numerical observation for Ly = 2000) and 
L ~ 2980 (theoretical upper bound from the mean field 
approach, Eq. H21|l . For practical numerical purposes 
at moderate system sizes, we note that such systems get 
into localizing regimes for x-system sizes L/t at one to 
two orders of magnitude than the previously derived up- 
per bound, Eq. if^ . 

Eventually, in the localizing regimes, we need to distin- 
guish between what was referred as diffuse or total local- 
ization in Section llVI total localization was defined as a 
case where the most probable break after departing from 
percolation, would happen in the zone referred as (1), i.e. 
corresponding to a distance r/i from the last burned fuse 
smaller than \J Ljl . Diffuse localization corresponds to 



cases where the next event would happen preferentially in 
zone (2), at moderate current perturbations, which corre- 
sponds to distances r jl from the last burned fuse ranging 
from \fLji up to a few Ljl. A criterion to distinguish 
numerically between these two regimes is thus to look, 
for a fixed large elongation Ly, whether the dependance 
of the average distance between successive events over 
the system width is such as d/^/Ll vanishes at large L, 
or on contrary remains finite or diverges. On fig. |H1 , the 
average distance between successive events was evaluated 
over 50 simulations and over the second half of the events 
before complete breakdown, which is in the localization 
regime for all cases probed. This average distance seems 
lowly sensitive to Ly = 200 or 2000 (20% difference or 
less between both sizes) , but the scaling as function of L 
shows that this distance saturates rapidly for a = 0.75, 
while it grows approximately as ^/L for a = 1.5. The 
extent over which this power-law corresponds to slightly 
more than a decade for L, which is the maximum achiev- 
able numerically since the elongation Ly has to exceed 
significantly L to be in the considered framework. This 
result is thus consistent with a transition from diffuse to 
total localization between a = 0.75 and a = 1.5 - the 
theory for Ly ^ oo predicts this transition at a = 1. To 
pinpoint more accurately the precise value of the transi- 
tion exponent (between 0.75 and 1.5 is not easy numer- 
ically: this would require a priori to look at the scaling 
of d{L) over more orders of magnitudes, for numerous 
a inbetween, which would represent a significant numer- 
ical cost and was not the main objective of this work, 
performed mainly as a numerical check of the analytical 
derivations carried out in the previous sections. 



VI. CONCLUSION 

The main results of the analytical calculation pre- 
sented in this paper are to be found in Fig. ^ There 
are three distinct phases of the fracture process depend- 
ing on the disorder exponent a and on the ratio between 
width of the lattice L and the lattice constant £ when the 
lattice is a cylinder of infinite length. The first regime 
is a percolation-like regime where the distance between 
successive failing fuses is completely random. In the sec- 
ond regime, named "diffuse localization" , the controlling 
parameter is L/t, while in the third regime, "complete 
localization" , the controlling length is t. 

To summarize the numerical results of Section 
regimes corresponding to percolation, diffuse and total 
localization have been clearly identified. The transition 
from diffuse to total localization is consistent with the 
predicted a = 1. The transition from percolation uni- 
versality class to localizing regimes is seen to happen un- 
der increase of either the exponent a of the quenched 
disorder power-law distribution, or the width of the sys- 
tem, as predicted by the theory. Nonetheless, the tran- 
sition width for a given exponent was found significantly 
smaller for the systems of finite elongation studied than 
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for the infinite elongation system considered analytically. 
This discrepancy was observed to be lower when the elon- 
gation Ly increased, and presumably corresponds to im- 
portant finite size effects. The transition from diffuse to 
total localization is consistent with the predicted a = 1, 
irrespectively of the system's elongation. 

Hence, in the diffuse localization phase, we expect a 
smooth variation of the damage profile scales L, while 
there is still no "strong localization" driving the burned 
fuse to merge at the lattice constant scale: region where 
arguments based on percolation in gradient should apply. 
In this case, the shape of the dipolar current perturbation 
(fig. lU, leads to most probable relative positions {x,y) 
of the next burned fuse relative to the last one, for which 
X and y are of the same order of magnitude, which is 
at the origin of the smooth quadratic maximum of the 
damage profile as function of y, as observed in This 
confirms the scaling arguments used to relate roughness 
exponent and correlation length divergence exponent in 
such systems. 



Appendix A: STATISTICAL LEMMAS. 

Consider p different types of random variables y, char- 
acterized by their cumulative distributions Pi{y) and 
probability density functions Pi{y) — dPi{y)/dy = Pl{y), 
for i — I. ..p. Next, consider an ensemble of ni ran- 
dom variables distributed according to pi, n2 accord- 
ing to '^p according to Up. In the limit where 
iV = X)r=i rii ^ 1, we wish to characterize 



(Al) 



The probability that some particular variable number j 
of type i, yij, would be equal to x, while all others are 
larger, is 



p,{x){i-p,{x)r-'l[{i-p,{x)r 



(A2) 



The probability that any of the variables of type i would 
be the smallest and equal to x, is the above times a factor 
rii. Thus, the wanted quantity may be written as 



/ xdxY^ 



M^){i-p,{x)r-'l[ii-p,{x)r 



cdx 



dx 



dx 



l[{i-pjix)r 



Setting p = 1 in Eq. ljA3|l . we have that 
dxil-Pi{x))^ . 



(A3) 



(A4) 



Since the function 1 — Pi ix) decreases continuously from 
1 to 0, for large n the product (1 — Pi{x))^ is equal to 
one for a; < Xc, up to a certain cutoff Xc, above which 
it becomes vanishingly small. The integral Eq. I|A4|1 is 
then simply equal to Xc- We determine Xc by invoking 
the standard saddle point approximation, which leads to 
the equation 



(A5) 



Pl{Xc) l-Pi{Xc) 

for Xc- By using rHopital's rule, p'l/pi ~ pi/Pi, this 
equation reduces to the condition NP{xc) = 1. Using 
m = Xr, we have Hal 



NPi{m) = 1 



(A6) 



Generalizing this result to p > 1 , we find by invoking the 
saddle point approximation for Eq. lABII , the equation 



St 



niPi(m) 



St 



-P,{m) 



(A7) 



If we now set 1 — Pj{m) « 1 as Pj{m) <C 1, and use 
I'Hopital's rule. 



ELi rijPj (m) J2^.^ njPj {r 



Eq. IjATp reduces to 



(A8) 



(A9) 



which generalizes Eq. (jA6|) . 

In the case of an infinite number of random variables 
indexed by a continuous parameter s, the equivalent of 
Eq. is 



n{s)P{m, s)ds — 1 



(AlO) 



The probability Qi that the minimum variable would be 
of type i (for a discrete set of random variables) , is from 
the above 



Qi = j dx 



n.,p,[x){l~P,{x)r-'\[{l-Pj[x)T 

(All) 

The same argument shows that (1 — Pi(a:))"* ^ Ilj/i (1 ~ 
Pj (x))"^ is equivalent to 1 when a; < m, or when x > m. 
Thus, 



= dxri^ptix) 

J x—~oo 

= UiP-Xm) . 



(A12) 

For a continuous set of random variables, the probability 
that the minimum variable would correspond to an index 
s £ [si, S2] is then 

n{s)P{m,s)ds . (A13) 
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Appendix B: DENSITY OF CELLS PER LEVEL 
OF SUSTAINED CURRENT 



To compute n(s), we will use an explicit expression 
Va{u) of the contours of iso-perturbation, shown in Fig. [21 
defined implicitly as 



f{u,Va{u)) = a. 



(Bl) 



Inverting this expression from Eqs. H10lljl(l . comes the 
following: if a < 0, there are two points of abcissa u 
satisfying this. Their ordinates are 



v^{u) = acosh 



2a — 1 , , , /cos^(u) sm^(u) 
cos(w) ± ' ' ' 



2a ' ' V 
These functions are defined for 



(B2) 



with = acos(4a/(4a - 1)). If < a < 0.5, there 

is a single defined function Va{u) satisfying Eq. I)B1|I 
for any u, which is the positive alternative of the above 
Eq. (|B2p . Last, if a > 0.5, the expression of Va{u) is 
identical, but once again this function has a finite support 
u mod [27r] e u^J*"^], with = acos(l - 1/a). 

Some examples of these auxiliary quantities were set in 
Figure |21 

To compute n(s), we reformulate the condition j G 
[je(l + s), je(l + s + As)] as f{u,v) € [a, o + Aa] where 
a = 2L'^s/{tt^F) and Aa = 2L^As/{Tr^l^) from Eq. ®, 
which through a Taylor expansion of f{u, v) in v around 
Va{u) defines 

n[s,As) [Va(u)-Wa{u)Aa,Va(u)]}{m) 



df 

Wa{u)= -l/ — {u,Va{u)) 

ov 



(B4) 



and thus n{s) 
as 



n{s) 
where g{a) 



^ I(x,y)Gn{s,As) ^^^^2^ '^ai^ expressed 
2L^ 



min(u 



g{2L^s/TT^f 



(B5) 

Wa{u)du. (B6) 



M = 



This number density g{a) was numerically evaluated from 
the above and displayed in Figure |2| We will also derive 
below the asymptotic behavior of this function around 
the special values a ^ 0, infinitely away from the burnt 
fuse, a +00, around for the near neighbors of the 
broken fuse, and a ~ 0.5, around the saddle point 
[u — TT,v — Q). As will be shown straightforwardly, these 
asymptots, displayed in Fig. |31 are: 



5(a) 
g{5a + 1/2) 

9{o) 



-a^o V2|a| (B7) 
^Sa-.o -21n(|5a|) (B8) 
^a^±oo l/2a^ (B9) 



Indeed, from Eq. itTUIl 

Wa{u) = 



(BIO) 



(cosh(wa(u)) — cos(u))'^ 



sinh(ua(it))[cos2(u) + cos(u) cosh(tia(M)) — 2] 



where Va{u) is given by Eq. (jB2|) . 

Developing the above around a '--^ 0"^ to main order 
in a is a direct exercise which leads to Wa(w) ~ —1/a 
when u g [7r/2,7r], and |wo(w)| <C 1/a for u g [0, 7r/2]. 
The integration of Wa{u) in Eq. (|B6p leads then to the 
asymptotic value in Eq. (|B7|I . 

Conversely, around a ~ +oo, f{u,v) ~ 2{u'^—v'^)/{u'^+ 
t)2)2 _ 2coa{26)/r^ where . Then, defining 



I{a) 



dudv. 



u>0,v>OJ{u,v)<a 



we have 



g{a) — —dl{a)/da 
and in polar coordinates 



I{a) = 



r/4 /-/ 



rdrde = I /2a 



(Bll) 



(B12) 



leading to the asymptotic form g{a) ^ \/2a}^ Eaa ljB9|) . 

Eventually, around 8a = a— 1/2 ~ 0, we refor- 
mulate the condition f{u,v) € [1/2 — Sa, 1/2] by a 
Taylor expansion of / to second order in v: Sf ~ 
f{u,v) - 1/2 = f{u,v) - f{u,V().5{u)) dyf{u,vo.5{u))- 
Sv + dvvf{u,VQ,c^{u))-6v'^/2. The first term dominates if 
Sa < \2{dvf)^/dvvf\ and then < Sf < 6a is equivalent 
to < Sv < -Sa/dyf. When Sa > 2{Syf)^/dyyf, we 
will have |(5/| < (5a if < < \/Sa/dyyf. From Eqs. 

IIB2IB10t . vo.5{u) = V2-cos2(u), and 5,/(w,wo.5H) 
has a single zero in m = tt (saddle point), dyf^ (tt — u)/4, 
whereas dyyf{u,VQ,z{u)) remains finite, tending towards 
a finite value a when u -k. Thus, we have 



/ 

J u,v 



dudv 



//(ti,t))e[0.5-<5a,0.5] 



Sa I VO{5a) 



^2(5aln((5a) + 0((5a) 



(B13) 



(B14) 



Differentiating this expression with respect to Sa shows 
that (7(0.5 — Sa) ~ — 21n((5a), for Sa 0+. Similar argu- 
ments lead to the same asymptotic form for a negative 
Sa, i.e. to the asymptotic form in Eq. ljB8|l . 



Appendix C: SPATIAL DIVISION OF THE 
INTEGRAL CHARACTERIZING EXTREME 
STATISTICS 

Here, we present the details of the subterms interven- 
ing in the implicit Eq. (jSJ, which determine both the 
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region of most probable next break with respect to the 
previous one, and the most probable jump in external 
current to reach this next break. 

We have to solve the implicit Eq. jSJ, with a derived 
distribution 

P{m,s)= \_>.^ ^' Ue{m{l + s)-j,). (CI) 

^ Je 

coming from the assumed quenched disorder distribution, 
P{t) — on < t < 1, with a > 0, transformed into 
a conditional probability of breaking for fuses having so 
far survived, using Eq. |0J). 

In the integrand of Eq. jSJl, 7i(s) is expressed as func- 
tion of g{a) through Eq. (|11|1 . and we approximate this 
last function by its asymptotic forms Eqs. ljB7IB8IB9|) 
around the singularities and in the tail of highest pertur- 
bations. 

We divide the support of the integral in seven zones: 

(1) : in the neighboring cells carrying maximum cur- 
rent, such asr/i < c where c is a finite large number. The 
nearest neighbors on the sides of the last broken cell carry 
the maximum current, Sj/je = f{2Tri/L,0)-n'^£^/2L'^ — 
1/4, corresponding to an upper cutoff um — L"^ / l"^ . 
Close to the origin, the current perturbation falls as 
l/r^, and this close zone is then defined by the condi- 
tion om/c^ < f{u,v) < um- 

(2) : The cells carrying a current such as 3/4 < 
f{u,v) < c^/27r^. On the x— axis, from Eq. ((Tn|l . these 
conditions correspond io2L/c < x < acos(--l/3)i/27r ~ 
L/3. We require that these two defined first zones have 
a common boundary, which by equating the cutoffs of / 
sets = Lji. 

(3) : A zone around the saddle point, defined as 1/4 < 
/Kt;) <3/4 

(4) : A zone such as Om < /("«,«) < 1/4. This zone in- 
cludes the far-range from the last broken fuse, on which 
the current is slightly increased by its presence. The 
lower cutoff is determined by the y— extent of the 
system (length of the band), which was so far omitted. 
When the finite aspect of Ly is taken into account, the 
current perturbation derived from boundary conditions 
at infinity, Eqs. (|9I1U|) is still valid in boundary condi- 
tions corresponding to setting the global current through 
the top and bottom plate, i.e. 

/ y-3{x,±Ly/2)dx/L^jo. 

Jx=-L/2 

Indeed, we find for Ly ^ i, that 

f{u,TTLy/L) ~ — cos(7t)/cosh(7r_Lj^/L) 
~ — cos(7i) exp(— 7rLj,/L) , 

whose integral is zero for u G [— tt, tt]. However, when 
counting the number of cells sustaining a given cur- 
rent, the condition \y\ < Ly/2 should be added to de- 
rive n{s). From the above, n{s) is unmodified when 
a > am = exp(— 7rLj,/L), or when a < —arm but in 



the neighborhood of zero for — < a < n{s) has to 
be modified according to 



g{a) = - (acos(-a/am) - tt/2) . (C2) 
a 

This function is pair, and increasing for a > from 
g{0) = l/a„ to g{a,n) = Tr/2a,n- 

(5) : A tiny zone of vanishing current perturbation, 
—flm < /("", v) < am, for which n(s) has just been deter- 
mined. 

(6) : A far-range zone where the current is screened by 
the last burnt fuse, —1 < f{u,v) < —am- 

(7) : A zone of highly screened current, —am < 
f{u,v) < -1. 

The zones corresponding to regions defined as flc, 
and 51/ in Section Hvl are respectively regions (1), (2-3) 
and (4-7). They correspond to regions, which are with 
respect to the last break, either close to it, at distances 
comparable with the lattice step, or "diffusively close", 
at distances comparable to the width of the system Ly, 
or "far" , at distances comparable with the system size 
Ly 3> Lx- 

Dividing here in seven subregions, will allow to use 
asymptotic forms of the seven corresponding subinte- 
grals. 

We will classify the regime of the system according to 
the zone where most of the integral in Eq. ^ is real- 
ized: We expect this zone to be either (4-6), in which 
case the system remains in a diffuse regime where there 
are no noticeable spatial correlations in the pattern of 
burnt fuse, or (1), which signifys the onset of a complete 
localization regime where the damage will develop in a 
concentrated zone scaling as the lattice size £, and tear 
through the system starting from this smallest scale, or 
(2-3) , which would denote the onset of a diffuse localiza- 
tion regime, where the characteristic distance d between 
the burnt fuse scales as L, the system's width. 

Noting e — 7r^£^/2L^ and A — m/je, where je,™ are 
the values of the external current at last break and for 
the next one, Eq. (O is then equivalent to 



EMA) = ^^^, (C3) 

where with 

p(a,A) = He[A(l -l-ea) - 1] , (C4) 
the /i-functions are defined from the above zones, asymp- 
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totic behaviors of n(s) and q.d. distribution as 
' (1 + ea)"A" - 1 



hi{\) 
/14(A) 



2a2 



-p(a, X)da, 



1/4 , 

21n( — )[(l + £(0.5 + <5a))"A^ 

-1/4 

p(0.5 + (5a, A)c?(5a, 
^Z"* 7r[(l + £a)"A° - 1] 



2|a| 



-p(a, A)(ia, 



(C5) 

l]x 
(C6) 
(C7) 



- (acos(— a/flm) — 7r/2) p(a, A)c?a (C8) 
a 



and /12 has the same integrand as /ii, but a support re- 
spectively on [3/4, c2/27r2], /ig has the same integrand as 
/i4 and bounds [—1, — a„j], /17 has the same integrand as 
h\ and has a support [— oaj , —1]. 

We first note that the term on the right hand-side of 
Eq. HU3|I is in the first stages of the process a large 
number, since the first fuse burns at an average value of 
the external current ji such as jf — 1/Nceiis = (-"^ /LLy, 
so that 2e(l — jf ~ 'K'^Ly/L > 1. It becomes a 
number of order e when ~ 1/2. 

We first look for solutions corresponding to an increase 
in external current, i.e. we look at the behavior of the 
/i— functions for A > 1- The Heaviside terms can there- 
fore be neglected for hi, ... , /15. Expanding the inte- 
grands to first order in e, and keeping the leading orders 
in c = ^JlJI., we obtain 



hi(X) - 


e[s{a)\" - 2n{c^ - 1)] 


(C9) 


s(a) = 




(CIO) 


/12(A) = 




(Cll) 


/13(A) ~ 


(ln(4) + 1)(A° - 1) + yA- 


(C12) 


/14(A) ~ 


-ln(l/a„)(A"-l) + — A" 


(C13) 




[|ln(l/a™)(A"--l) + ^(l- 


i)A"] X 




He[A(l + ea„0-l] 


(C14) 



The other /i— functions can be shown to be negligeable in 
front of these, and are not displayed here. In the early 
stages of the process, the first terms of /14 and will 
dominate, i.e. the singularity of n around s ~ corre- 
sponding to the furthest zone of the last break will be 
preponderant, and the threshold to next break will be 
set by 



7rln(- 



1 



1) 



)(A^ 
i.e. A"-l = 
which corresponds to Eq. (|15|> . 



2e- 
1 



1 



Nr.. 



dlsJe 



(C15) 
(C16) 



This leads to a second break happening on average 
when j2 = Aiji = (1 + l)i/"ji = {2/N,M,y/''. As long 
as the process remains in a diffuse regime where burned 
fuse are far away from each other, at distances compara- 
ble with the y— extent of the system Ly, we can go on 
with this mean field theory, and use the same arguments 
to evaluate the probability of finding the next (n+1 — th) 
burned fuses with respect to the previously burned ones. 
In this situation, the fuses of interest are far from the 
already burned ones (/14 and dominate the integral), 
and thus have thresholds above je, the external current 
reached so far. We can then use the same arguments, 
replacing the average current at first break ji by the av- 
erage external current level at n— th break j„, to obtain 
by reccurence that j„ = {n/NceiisY^°' , which corresponds 
to Eq. H17|) . Indeed, by definition, jn+i — ^njn, and from 

Eq. (IHTeli . 



A^ - 1 = l/{Nr,Usj:) = 1/r 



(C17) 



Thus, A„ = [(n+l)/n]i/" and j„+i = [(n-f l)/7Veei/«]'/", 
establishing the result by recurrence. Note that this 
result is a simple consistency check of the present the- 
ory: In this percolation regime, the concentration of the 
current around the broken fuses is negligeable, and the 
level jn+i of (n -I- l)-th break should be given by the av- 
erage minimum of the thresholds of the entire system, 
under the condition that all of them have survived up 
to the current jVi- This formulates as P{jn+i\jn+i > 

jn) ^ Un +l - J")/(1 " Jn) = ^/Ncells, which is CXactly 

Eq. iHlBl) . 

A carcfuU analysis of the above functions shows that 
in the limit y, these first terms of /14 and /ig always domi- 
nates: they are larger by a factor ln(l/am) = Ly/L than 
the others terms proportional to (A" — 1), and all of 
the corrections proportional to e vanish. Thus, in this 
limit, Hf dominates Eq. H14|l . and no localization hap- 
pens. This lasts during the whole process (as long as 
the interactions between already burned fuses are weak 
enough for this mean field theory to apply), and in this 
limit of large disorder a — > 0, the process remains diffuse, 
in a percolation-like regime, up to the moment where 
— 1/2, i.e. when P{je) = 1/2, which corresponds to 
the critical percolation threshold. 

On contrary, for a —^ +00, s{a) is diverging faster than 
all prefactors of A" in the above, and /ii(A) defined in 
Eq. (|C9|) . i.e. dominates the Left Hand Side in Eq. HC3p . 
i.e. Hc{X) dominates in Eq. (|14|l . thus leading to a level of 
next break set by /ii(A) = 2e(l - j^)/j^, i.e. Eq.((THll. In 
this case, the next fuses to burn will be the ones carrying 
the highest current perturbation, i.e. the near neighbors 
on the sides of the first one (moreover, s(a) is dominated 
by the contribution between l/47r and 1/477/2^, i.e. in 
this limit a —i- 00, the nearest neighbors will burn with 
certainty) . 

In the range of finite a, the system can be driven to 
a third regime if /12 or dominate in Eq. (|C3|) : cor- 
relations in the damage start to be significant, but the 
characteristic distance to the preceding burnt fuses is in 
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a range between \fLl and L, and does not scale as the 
lattice constant t. this is the regime which we refer to as 
"diffuse localization" . 

We have shown above that in the early stages of the 
process, for a finite disorder, i.e. at finite a, the system 
starts in a percolation-like regime, i.e. that i/c (h^ and 
h-i) dominate in Eq. (|C3|I . This can be the case up to 
the percolation transition, when Pije) = 1/2, for large 
enough disorder, i.e. small a. But as the weakest bonds 
are broken, the process can transit to one of the two other 
regimes: 

either the "total localization regime", where rupture 
proceeds via jumps between successive burned fuse, 
whose size d is close to the lattice spacing, d ^ 
. Either the "diffuse localization regime" , characterized 
by a distance d between successive bonds scaling typ- 
ically such as \f{L€) <^ d <^ Ly when L ^ co and 
Ly ~t oo. This transition towards the diffuse localiza- 
tion regime should happen when H^i stops being small 
compared to He, while Hf remains negligeable compared 
to these. Technically, this transition can be determined 
by equating Hc{\) — H(i{X), which corresponds to lead- 
ing order to /12(A) = ^4 (A) -I- /i6(A) (it can be shown 
that the contribution of the saddle point hs (A) is always 
negligeable in front of /12(A) as long as = L/£ ^ 1). 
From Eqs. I|C11IC13IC14|I corresponds to leading order 
in 1/Nceiis to the condition 

2M^)=2-|^' (C18) 



which corresponds to Eq (|20|l . As detailed in Section ITVI 
if this condition is not met at the percolation transition 
when j'g ~ 1/2, the system stays in the percolation like 
regime, whereas if this condition is met before and He is 
still negligeable, the system transits towards the diffuse 
localization regime before breakdown: from Eq. (|('18p . 
this transition never happens if a < am = i^TfXyi)' which 
corresponds to Eq. (|^ . 

To characterize the boundary between transition to to- 
tal localization, and diffuse localization, i.e. to obtain 
Eq. (|23|l . we have considered whether indeed Hc{X) is 
still small compared to Hd{\), with A approximated by 
its value at percolation threshold, i.e. fixed by Eq. ifT^ 
with = 1/2. Expressing in this way the condition 
hi[X) = /12(A), with the detailed expressions in Eqs IjCQI - 
ICllfl . together with the definitions am = exp(— ttLj^/L), 
e = 7r2^V2i^ and = L/i, leads to Eq. 

The estimate of the current value, Eq. H24|) . where the 
system transits towards total localization when a > 1, 
is obtained by equating /ii(A) = /14(A) -I- /i6(A), and by 
noting that the last term of /14(A) dominates the right 
hand side. 

Eventually, Eq. H25() . which determines the boundary 
between the percolation like regime and the total local- 
ization regime in the a-L/t space, is obtained by equat- 
ing /ii(A) = /14(A) -|- /i6(A), with A evaluated from its 
value at percolation threshold, i.e. fixed by Eq. l(T5)l . 
with j"^ = 1/2. 
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